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(54) 

(57) The signal to be measured is converted into 
digital form and read into a memory 13 using a tngger 
as the reference, and a signal of a length 1 15 that of the 
signal to be measured, in the vicinity of the trigger, is 
Orthogonal-transformed (14). A frequency erra co, of 
the transformed output h.Qi is estimated (15) and is 
used to correct the output l,.Q, (16). The corrected out- 
put I2.Q2 is subjected to processing by a receiving filter 
(17) and the filtered output Is.Qa is synchronized with a 
PN to obtain an ideal signal R,.Ro and a phase dif- 
ference XI is calculated. The signal I3.Q3 and the ideal 
signal R,.Rq are used to obtain a frequency error ^ 
and a phase error ^ of the signal b-Qa- AH signals to be 
measured, stored in the memory 13, are subjected to 
orthogonal transformation (23) so that the frequency 
error tDi+oog and the phase error ti+t2 are removal, 
and the transformed output Ig.Qs and the signal R,.Rq 
are used to calculate modulation parameters as is the 
case with the prior art. 
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Description 

BACKGROUND OF THE INVENTION 

5 The present invention relates to a method for measuring modulation parameters of the transmitted signal quadra- 

ture-modulated by a preknown code (data), such as a carrier frequency error, the modulation accuracy, an IQ origin off- 
set, an amplitude error, a phase error, the waveform quality artd a time alignment error, as is the case with the 
transmitted signal from a base station defined by, for example, the standard of CDMA (Code division Multiple Access) 
digital mobile communications (EIA/TIA/IS-95 Mobile Station -Base Station Compatibility Standard for Dual-Mode 

10 Widetjand Spread Spectrum Cellular System, which will hereinafter be referred to as IS-95 standard). 

The output signal from a conventional transmitter of the QPSK modulation system needs to be demodulated accu- 
rately on a bitwise basis. On this account, criterion for assessment by the standard are rigid ^nd the output signal from 
the actual transmitter to be measured is also a relatively clean signal of a waveform not seriously distorted or disturbed, 
that is, an easy-to-demodulate signal. Since the signal of the system of the IS-95 standard is a spread spectrum signal, 

15 however, it is not required to be demodulated accurately bit by bit. That is. the QPSK signal presents no problem even 
if it is a signal of a distorted or disturbed waveform such that a bit error arises when it is demodulated after quadrature 
detection. Hence, it is necessary that a parameter measuring method for determining the waveform quality of the trans- 
mitted signal from the base station, which meets the IS-95 standard specifications, be able to determine the waveform 
quality even if the waveform is distorted or disturbed more than in the case of the QPSK modulated signal. 

20 In view of the above, p is defined by the IS-95 standard as a waveform quality parameter, not as the modulation 

accuracy. The parameter p is obtained by normalizing the value of correlation between the modulated high-frequency 
signal from the transmitter and an ideal modulated signal by their powers. When the both signal match, the parameter 
p assumes the maximum value 1. 

In the CDMA system the coincidence of timing among base stations in the communication system is indispensable 

25 because of the property of the spread spectrum signal. As regards the transmitted signal from the base station, it is nec- 
essary that the timing for sending the high-frequency signal coincide with certain timing of the standard time in that base 
station. The tolerance of this timing is also standardized by the IS-95 standard. This timing is defined by the difference 
between a trigger that is sent from the base station and spedfic timing of a PN pattern of the high-frequency signal 
transmitted from that base station. This trigger is sent at the beginning of the PN pattern period in the standard time in 

30 the base station. This timing error is called a time alignment error 

The CDMA signal is spread by a high spreading chip rate rather than by the transmission bit rate. Hence, in the 
assessment of a modulated signal of a duration corresponding to transmission bits of the same number as that of the 
bits subject to assessment in the conventional standard for measuring the modulation accuracy of the QPSK modulated 
signal, the number of chips will be appreciably large even if the duration is the same. In the conventional scheme for 

35 measuring the modulation accuracy of the QPSK modulated signal using the chip number as the symbol number, the 
number of symbols to be measured is large and the corrputing time is long accordingly. By the way, in RDRSTD-28 that 
is the PHS standard in our country, a signal of a 1-burst-600-jiS period is subject to assessment and the symbol rate in 
this case is around 200 kHz, whereas in IS-95 the chip rate is approximately 1 .2 MHz. If a signal of the 600-fiS period 
is subject to assessment as in the case of the PHS standard, the number of symbols to be estimated will become about 

40 six-fold. 

In the measurement of the waveform quality, it is standardized that the transmitted signal to be measured is applied 
to a complementary filter to eliminate intersymtx)! interference that is caused by a filter and a transmitting phase equal- 
izer at the transmitting side. 

To this end, the complementary filter is defined to provide a Nyquist filter characteristic when the sending filter, the 
45 transmitting phase equalizer and the complementary filter are connected in cascade (connected in series). 

In the CDMA system the use of such a complementary filter characteristic in a filter at the receiving side provides 
the optimum receiving characteristic. Consequently, this complementary filter serves as a receiving filter of the meas- 
uring apparatus. 

The complementary filter serves as a narrow band-pass filter for the high-frequency signal but a low-pass filter for 
so the base band signal. In the measurement of the waveform quality, processing by this filter is carried out in accordance 
with the carrier frequency of the input signal, but in this instance, if the carrier frequency is not accurate, there is a fear 
of the input signal being partly cut out because the filter is narrow-band. 

Further/the receiving filter is equippefdwithr group delay frequency characteristics that Cancel those of the transmit-'' 
ted signal. On this account, it is necessary in the measurement of the waveform quality, too. to measure the timing 
55 (phase) of the Input signal after the filter processing. In this case, a major or serious error arises in the timing measure- 
ment unless the carrier frequency is accurately estimated. 

Conventionally, this filter processing is repeated a number of times to repeat the estimation of the carrier frequency 
as disclosed in the specification of U.S. Patent No. 5,187.719 (issued february 16. 1993), especially in Fig. 15 and the 
associated disclosure. The timing (phase) of the input signal is estimated only once prior to the filter processing since 
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"^r:rer=:~i.^— ^^^^^ 

spreading code of the measuring apparatus 7*^^*^^.^^^^^^ standard is called a pilot channel. The 
one Of Channels for the t^«"--'«^!|"^' ^^ul^^^^^ ^'^ 

'^""^^s^SKS^^^^^ 

ternary in the prior art to convert or transforrr, the '"P^J J^^^^^^^^^^^^^ Se (PN c^ef of the pilot signal, i.e. a quad- 
trequency set eighttimes ^^'9^-*^^^^^^ the correla«on between them while 

ruple oversampling. Since the spreading , '° ^ operation is appreciably time-consuming. 

-rtrm=ro:e=raSE^^^^^ 

S^pr^bleband Signal. This willbedesc^^^^^^^^^ 

For example, a QPSK complex base band signa « '!^'"P'^^^^^^^^ ^ by the oscillation output from a volt- 
and a differentiation circuit 11 3. In the sampling circuit 1 ^^s'J"^^^^^^^^^ a,so sampled by the out- 
age-controlled clock (VCC) generator 1 1 4 and the output^om Je ^ ^ tHe sampling 
p'ut from the vo«age-contro.led clock 9^"^^ ^ ^^^^^e w^e n ^^^^^^^^ ^yTn ideal signal (a reference sig- 
circuits 1 12 and 120 are fed to multipliers 1 ^ ^ and 1 m wherem ttiey a P ^ ^^^.^ 
nal) exp(-iek) (©^ being a k-th phase of dernodj^ate^^^^ ,,9^ namely, they are 

with the symbol timing. ^^in^iion of the svmbol timing of the QPSK complex base band signal by 

them\^rmlS:itLS^^^^ 

'"t^'^Sac!: a logarithmic likelihood function A, can be set as follows: 

A = Re[exp(-j4.)fr(t)R*(t-T)dt] 

s,^ro=:^N"=rsii" "s- ^ ^ — 

From the above equations x is calculated to satisfy the following equation. 
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(1) 

Re[Z(T) • (a2*(x)/ax)] = 0 

where 

Z(t) = fr{t)R*(t-T)dt 
where I indicates an integration from 0 to T and 

(3) 

R*(t-T) = i:g(t-T;-kT^-)exp(-jek) 
spending to T Substituting them into Eq. (1) gives 
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Re[Z(T) • dZ*(i)/ch] = Re[5:{exp(-jek) • Y ,,(t)} • Z{exp(-j0k)c?Y ^,(T}/dx]] = 0 (4) 

where 

5 Y,(t) = Jr(t)g(t-T-kTJdt (5) 

In the above. I is an integration from 0 to T 

That is, 5;{exp(-jek) • Yj^(t)} corresponds to the output from the circuit 118 in Fig. 1 and 2:{exp(-jek)c)Yk(x)/^} corre- 
sponds to the output from the circuit 119. The VCC 1 14 is controlled so that the real part of the products of these out- 
10 puts is reduced to zero, and the sampling timing of each of the sampling circuits 1 12 and 120 is estimated from the 
symbol timing of the input QPSK complex base band signal. 

The configuration of Fig. 1 is to estimate the symbol timing by the maximum likelihood estimation method. 
The Fig. 1 configuration is formed by hardware, and in the case of implementing it by software, no highly accurate 
solution could be obtained because the input signal is a discrete time signal. To obtain the solution with high accuracy, 
15 the sampling rate (the sampling speed) must be increased -this inevitably increases the amount of processing by an 
interpolation filter and hence increases the processing time. 

In the estimation of the carrier frequency in the conventional parameter measurements, in the case of a four-phase 
(n=4) PSK signal, that is, in the case of the QPSK signal, the intermediate-frequency (or high-frequency) QPSK signal 
in the digital form is raised to the fourth power to thereby remove the modulated signal. That is, the QPSK signal 
20 assumes any one of four phase points sequentially displaced ti/2 apart in accordance with the modulated signal, so that 
by raising the QPSK signal to the fourth power, any phase states become integral multiples of 2n and the modulated 
signal is removed. 

The signal thus raised to the fourth power is subjected to fast Fourier transform, then a frequency component that 
provides the maximum peak is extracted from the result of the Fourier transform and frequency divided down to 1/4, and 
25 the frequency of this divided output is estimated as the carrier frequency of the input QPSK signal. This estimated fre- 
quency is used to convert or transform the input QPSK signal to a complex base band signal through Hilbert transfor- 
mation or orthogonal detection, besides the estimated frequency is used to estimate the deviation of the QPSK signal 
from the standard carrier frequency. 

The input QPSK signal is a digital signal, and if its sampling frequency fg does not satisfy a condition 2fm<fs with 
30 respect to the maximum frequency fm of the QPSK signal, aliasing occurs and original information cannot be retained. 

Since the QPSK signal is raised to the fourth power, its maximum frequency is four-fold, that is, 4fn,. If the sampling 
theorem is not satisfied in this instance, aliasing occurs and the carrier frequency cannot accurately be estimated. 

It is therefore an object of the present invention to provide a method with which it is possible to correctly measure 
modulation parameters of the input signal evoi if it is prone to produce a demodulation error. 
35 Another object of the present invention is to provide a method which permits accurate measurement of modulation 
parameters of the input signal with a small computational complexity and in a short time. 

Another object of the present invention is to provide a method which enables modulation parameters of the input 
signal to be measured by a signal of a shorter duration. 

Another object of the present invention is to provide a modulation parameter measuring method which permits 
40 rapid synchronization of a spreading code with the input signal. 

Another object of the present invention is to provide a modulation parameter measuring method which enables 
symbol timing to be estimated by software in a short time. 

Still another object of the present invention is to provide a modulation parameter measuring method which permits 
an accurate estimation of the carrier frequency. 

45 

SUMMARY OF THE INVENTION 

According to the present invention, a input digital quadrature-modulated signal is subjected to orthogonal transfor- 
mation to a first complex base band signal In a first step. In a second step a frequency error, an initial phase error and 

so a timing error of the first complex base band signal are estimated and an ideal reference signal is calculated which is 
synchronized with the first complex base band signal. In a third step the input digital quadrature-modulated signal is 
subjected to orthogonal transformation to a second complex base band signal so that the frequency error, the initial 
phase error and the timing error are corrected. In a fourth step a frequency^error^and an initial phase enror of the second' 
complex base band signal is estimated from the second complex base band signal and the ideal reference signal. In a 

55 fifth step the estimated frequency error and the initial phase error of the second complex base band signal is corrected 
to provide a third complex base band signal. In a sixth step a modulation parameter is estimated from the third complex 
base band signal and the ideal reference signal. 

In a seventh step the frequency error of the first complex base band signal estimated in the second step is cor- 
rected to provide a fourth complex base band signal, which is subjected to processing by a receiving filter or comple- 
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.i*7.*s I «al ;;:S™e"rn5S „^ >h, « convex 

W » second step, symbol decs,™ ports o^^ ^aaed Iron, » M cowta 

band signal seqosocs, then samples nea, the estirrated ^T^l^L ^SlS sheading oodes is calculated. In 

base band signal s«,uenoe, and the ~"«»««^ l"'*^ 1^ ''^^'=^2,°^^^ o, imaginary part 

this instance, the correlation calculation .s conducted only lor the real parnin pnas 

used to calculate the symbol timing. 5i , ^915^^ to the n-th 

in the second step, the first complex base band ^'g-^j' '^ra^^ to the n th P ^ ^^^^ ^^3^^ 

power is subjected to discrete Fourier transform, and '"^^"^"'^^ .^^^^^^ 3^ tt,e divided frequency and the 
Lsform. The frequency of the --J^ P^j^^^^^^^^^^^ o ei-te the carrier frequency of the 

ZX^^^^^ ;2S^thT.rra; sati.ie^s n=.^ (where m.. a. ...)• 

RRIFF DESCPiPTiniM OFTHF DRAWING S 

Fig. 1 isablocKdiagram showing thefunCona, ^-'^^^^^^-'^^^^ method 
Fig. 2 is a block diagram illustrating the functional ^^^''S" f the present invention; 

formeasuring modulation parameters of a digrtal ^"^*f " ^^"'^^^^^^^^^ ^5 in Fig. 2; 

35 nFRr-RiPTION nP thf PREFERRFD FMB ODIMENT 

l„P,.a«re.illus.r«ed^.^2^c».^----^^ 

X-'a^-^rpSaTe'SroKp^^^^^ is call. a1o« »de. and 
rc™^^be«een.He»»™pa~^ 

The pilot signal (hereinafter referred to as a BS P"o* ^'9"^''/'^*"^'"'". .^rm\nB\ 1 1 into an AD converter 12, 
tion is co'nverteS to an intermediate-frequency signal, ^j^";^ J^^^^^^^^^^^^ ote stored in a memory 

wherein it is converted to a digital signal, and the digits s gnal o t^e P^^'^^^^^^ ^^^^ ,,^„ency of a spreading 
1 3. The sampling frequency in the AD converter j 2 is set '° ^« hat whenTwrite or read reaches 
PN pattern of the BS pilot signal. The memory 13 '^^^^'^'^J.^^^^^^^^ J^ress increment and write halting 

the last address, it starts again at the address 0 and •'e write in rnemo^ an .^^^ ^^^^^^ 

are controlled by a trigger signal that is fed from the o.^'de^ ThaUs_ the ^^9.^1 s.g"a,s ^^^^^ 
1 3 are befae and after the trigger signaL Some ^ ^/Xu, before and aSer the trigger signal. 

2048 samples (in terms of the chip number. 1/5 i.e. about ^^0° ^^^r^^tlt smaHer than the data length of the period for 
in this instance, the number of samples of the digita^ l"!"'^jSn^^^^^ but it is true that the dat^- 

measur«v«,^ti^ined by the ^-^^^'^^^^^''^''Zt^Zs^^ to an orthogonal 

rsrei^p^rst^^^^^^^^^^ 
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complex base band signal is raised to the fourth power in a fourth power calculating section 124 to degenerate the mod- 
ulated component to obtain one phase state, then the signal raised to the fourth power is subjected to fast Fourier trans- 
form by an FFT section 125, and that frequency in the FFT output at which the maximum peak is obtained is frequency 
divided down to 1/4 in a frequency dividing section 126 and the standard frequency of the CDMA digital cellular tele- 

5 phone system is added to the frequency-divided output in an adder section 127 to obtain the carrier frequency. The 
sampling rate of the first complex base band signal is 4.9152 MHz, four times higher than the chip rate, the number of 
FFT points in the FFT section 125 is 1024, frequency resolution is 4.8 kHz and the frequency estimation accuracy is 
±2.4 kHz. The output frequency of the frequency dividing section 126 represents a deviation from the carrier frequency 
defined by the standard, i.e. an estimated frequency error. Since the data used for the estimation of the frequency error 

10 is data neightx)ring the beginning of the PN period, even rf it is displaced from the period to be actually measured, the 
carrier frequency or timing error does not change in a short time in the transmitter, and hence it is regarded as the car- 
rier frequency error in the measuring period. 

The thus estimated frequency error <o^ is used to correct the first complex base band signal \^,0^ in a frequency 
correcting section 16 (Fig. 2). This correction is made by the following calculation. 

15 

'2 +jQ2 +jQi)exp(-ja)it) 

The thus corrected complex base band signal l2,Q2 is subjected to delay complementary filtering in a receiving filter 
17. That is. In the base station the signal to be transmitted therefrom is subjected to filter processing so that it has a 
20 characteristic which will cancel the group delay frequency characteristics in the mobile station to which the signal is to 
be transmitted. For this reason, the receiving filter 17 cancels the group delay frequency characteristics of the filter in 
the base station. 

The filtered complex base band signal I3.Q3 is provided to a PN synchronized ideal signal generator section 18. 
wherein it is synchronized with a locally generated PN pattern. Prior to the synchronization, symbol decision points are 

25 estimated. The synchronization of the entire sampling data (in this example, the number of samples is eight times as 
large as that of one chip) vwth the PN pattern data inevitably increases the computing time. To avoid this, in-phase and 
quadrature PN codes are used to pre-estimate symtx)l decision points of the QPSK modulation and only sample points 
closest to the symbol decision points, that is. data every eight samples, are synchronized with the PN pattern data. The 
symbol decision points are estimated as described below. As shown in Fig. 4. the in-phase and quadrature components 

30 of the corrected complex base band signal l3(k),Q3(k) Is squared In multiplier sections 132 and 133, respectively, and 
the squared outputs are added together in an adder section 134 to obtain an annplitude squared sequence of the com- 
plex base band signal l3(k), Q3(k). which Is subjected to discrete Fourier transform in a discrete Fourier transform sec- 
tion 135. A symbol (chip) frequency component in the transformed output Is extracted and the arctangent of the 
extracted component is calculated in an arctangent section 136. From the arctangent output a sequence of zeroes is 

35 extracted, and the zeroes are alternately separated into two sequences, the one of which is estimated as a sequence 
of symbol decision points. 

In the above example the sampling frequency of the complex base band signal sequence is chosen to be eight 
times higher than the chip frequency In this instance, a sequence of i3(k)^-i-Q3(k)^ is multiplied by eight cosine values, 
for exarrple, 1, 0.707. 0, -0.707, -1. -0.707, 0 and 0.707 in a sequential order and a sequence of li(k)^+Qi(k)^ is simi- 

40 larly multiplied by eight sine values one after another. By this, the chip frequency component is obtained which is the 
result of the discrete Fourier transform. 

The samples near the chip decision points estimated in a symbol decision estimating section 131 are decimated in 
a decimating section 137 from the real part (the IniDhase component) l3(k) in the complex base band signal sequence. 
Assuming that samples in the sequence of the real part l3(k) are at the positions indicated by crosses X in Fig. 5 A and 

45 the estimated symbol decision points at the positions indicated by white circles Q. one of the samples at the posi- 
tions X which are closest to each position Q is taken oat as shown in Fig. 5B. In the case where the sampling fre- 
quency Is chosen to be an integral multiple of the symbol frequency, once the difference between a first symbol decision 
point and the sample point corresponding thereto is determined. It is enough to take out the subsequent samples at 
intervals of 4Ts (where Tg is the sampling period). The correlations between a sample sequence l3(4k) of the real part 

50 \^{k) near the first symbol decision point and spreading codes R|(4k) and RQ(4k) from a spreading code generating sec- 
tion 140 are calculated in correlation calculating sections 138 and 139. The results of the calculations are squared in 
multiplier sections 141 and 142, and the squared outputs are added together in an adder section 143. As the result of 
this, a correlation value A'(4k) bietwe6n thereal part l3(k) and the spreading codes R|(4k) and RQ(4k) is proCrided'from 
the adder section 143. The phases of the spreading codes that are generated in the spreading code generating section 

55 140 are shifted chip by chip so that the correlation value A'(4k) becomes maximum. When the correlation value A*(4k) 
becomes maximum, the spreading codes R|(4k) and Rq (4k) are synchronized with the input l3(4k). Incidentally, the cal- 
culation of the correlation value A'{4k) needs only to be conducted over ±1 00 chips of the real part l3(k). 

In the above-mentioned pilot signal, the correlation between the in-phase and quadrature PN codes is uniformly 
low, so that even if the phase of the complex base band signal is not correct, correct synchronization could be accom- 
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B,(4k) and RQ(4k) may also be similarly '^^^^^^^^^^ «1 a spread spectrum signal obtained 

m a different case, for example, in the case of synchron^ng PN c^^s wrt P ^ correlations 

by spreading the QPSK modulated ^'9-. by a ^prea^^^^^ ^^^^^ at 
between l3(4k). Q3(4k) and ^iW; RqW J" m^^^^^^^ 

intervals of the symbol sequence (41^ m the f "^j' ^^^.^^ This estimation can be made with accuracy 

is smaller than in the case of ^^Iculat-ng at .n er.al^^^^^^^^^ q ,he symbol decision 

^^Trp^yr™ 

Turning back to Fig. 2. the ideal signal R,.Rq thus P-^^^"^ "^^^^^^ ,3,q3 is compared in a frequency 

1 8 and the measured signal (complex base band ^'f^^^j^'^^J^JJ^^^^^^^^^^^ phase of the input signal and 

and phase error estimating section 19 to estimate a earner ' ^ corresponc«ng symbol decision points is 

me phase difference between the ideal ---1!^.^^^^^^^^^ That k representing time on 

also calculated. If the phase difference ' ' " is^ linear function of time and its intercept 

the abscissa and the phase diHerence on ^l^^ the inclination to the time axis is the earner 

or segment, that is. the phase difference a « "J.^^^^^^^^^^ by x and y. respectively, the incl.- 



b = (Sy , • XX i - nxx j • y |)/(2x , • Sx | - nSx , ^) 



35 



above equation becomes as follows. 

b = (5:y j • Si - nXi • y jVC^i • si - nXi^) 
where S indicates the summation from i=1 to n. Si l''}^^^^^^^^^^ nosrtillilJiquelc^^ error exists when the phase 

nal TiC" s foll^iSg equation in a frequency and phase correcting sect»n 21. 

I4 ■tp< •Os*P3)«>''<l»2t* e„) 

„ext*ean»t«.«d.a,^..n.,ses«.t^«a.2a^a.,J^^^^ 

rsn?roSrrs.^^trn:rsr;rn,Ca.on.,a— . 

A(,|,, f, x) = Const. {exp(-j4.)C{f, t) + expG<t>)C*(f, x)} 

(2) 

C(f, t) = 0 Z(t)exp(-2jiift)R*(t-x)dt 

• ^-,i,.,HT ie the measuring time for estimating the parame- 
so where Z(t) is the measured signal, R(t) is the reference signal and T„ .s the measur g 

*""The continuous signal R(t) and the discrete signal R(k) bear the following relationship. 

R(t) = R(kT3) 

where T, is the sampling period. disclosed, for example, in Proaks, "Digital Commu- 

The principle of estimation similar to that by ^^^^ ' ^"^^^ soecifically. the following simultaneous equations, 
nication." 2nd ed. pp. 333. Eqs. (4. 5 71). McGrow-^^^^ ^^.^^ p,,3^eters f and x. 

in which equations differentiated by the respective parameters are set ai 
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exp(-j<t»)C(f. t) - expG(t>)C*(f. x) = 0 (3) 

Bxp{'\^)dOldx + exp(j(t>)aC7aT = 0 (4) 

exp(-j(t>)aC/£)f + expG(t>)aCVaf = 0 (5) 

By simultaneously solving these equations, they can be modified to the following simultaneous equations without ^. 

d\0{l x)| ^/cH = 0 (6) 

5|C(f. T)|^/af = 0 (7) 

From Eq. (6) x is calculated setting f=0. Concrete steps of the above calculations will be described below. 
In the first place. C(0, t) is expressed with the integration replaced by the summation as follows: 

C{0. X) = T32Z(kT3) - R^(kT3-x) (8) 

where Z indicates the summation from k=0 to K-1 and Tg^T^j/e, Tc being the chip period. The measuring time 
Tq=KT3 . The continuous signal R{t) is expressed by the phase Op of the ninth chip as follows: 

R(t) = SexpOe p) • s(t - nTc) (9) 



where r indicates the summation from n=-oo to oo and s(t) is the characteristic of the base band filter In practice, the 
summation of Eq. (9) must be effected in a finite or limited period. To this end, the base band filter is assumed to be 
25 zero-phase and to have a response time Tp That is. 

s(t) = 0(|t|>TF/2) (10) 

In this case, the ideal signal can be expressed as follows: 

30 

R(kT^ - x) = S^expOe Js({8n-k}T3 + a) (11) 

Then, using Eq. (11), Eq. (8) can be expressed by the following equation by the use of the discrete signals Z(k) and 
R(k). 

35 

M 

C(0.x)=T3 X 2^2({8n-m}T3)exp(-jen)s(mT3+x) (12) 

nn=-M n 

40 Letting T p=(2M+l )T ^ , and the measuring time be represented by KTs (k=0 to K-1 ), the Op must be measured for a time 
n=- M/8 to (K+M-1)/8 . Next, s m(T^)=s(mT5+x) is approximated by a quadratic expression of x as follows: 

s^(x)=a^+bp,T + c^x^ (13) 

45 However, since t to be estimated is within the range of |x|<Tg/2 , the approximation needs only to be accomplished 
within this range. Using this, C(0, x) is given by the following equation. 

M 

C(0.t)=T3 2 2;Z({8n-m)T3)exp(-je JCa^+b^T+c^T^) (14) 

50 m=-M n 

^T3(A+Bx-i»-Cx^)" ' ^' ' (15)^ 

55 A. B and C in the above are given by the following equations. 

A=EEZ({8n-m}T Jexp(-je Ja ^ (16) 
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B=^;^Z({8n-m)T s)exp(-je „)b „ 

m n 

C=;^;^Z({8n-m}T 3)exp(-ie ^ 



(17) 



(18) 



Substitution 
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of Eq. (15) into Eq. (6) gives the following equation for x. 

Re[C(0. .)5C*(0, x)/^] = (T ^ Rel(A . Bx . Cx =^)(B- . 2C-x)l =0 (^9) 

™s.sacu.iceq.a«onfo.x,..ifar,ne.app— 
calculating the clock phase x is given as follows. 



15 ^ = -Re[ABn/{|B| ^ + 2RelAC*]) 



(20) 



This x is the aforementioned xg. het««een the ideal reference signal R,(8k)+jRQ(8k) and the rneas- 

That is. as shown in Fig. 6. the cross correlation ,51 . then the calculated result .s filtered 

ured Signal U(kHQ4(k) is calculated -^l-^^VT^Z rZ^^^^^^ 

by filters 152. 153 and 154 of filter ''-.f"f ""^^^^^^^ section 155 to calculate Eq. (20) to obtain 

and (1 8) to obtain A. B and C. which ^^^'"P"* ^j'^"^,^'^^^^^^^^^ for example, as described below. To prevent 
the tming error xj. The filter coefficients '^"j.^^^^'ll^ro ntSrf Ts impulse response to the first zero point needs 
SteUmSol interference in the Nyquist ""e^ ^-J^^^^^^^ S) feCa rSJ-off coe^icient be represented by a. the 
25 to be T, (the symbol period, in this example, the chip peno d).LeTi g 
impulse response of this filter is given as follows: 



sin(7 it/Te) cos(jtat/Tj) 
SW= ■ (ntn-,)"' 1 -(2atn-c 



30 



35 



40 



45 



50 



This impulse response is approximated by (13) ~ 

in an orthogonal transform, ^^^^"^"'^y' P'^^f " XiH^ the simal delay x. are used to perform orthogonal trans- 
so far. that is. the carrier frequency error, the '"Jj' P^J"^^^^^^^^^ stored in the memory 13, by which is 

formation of all the intermediate-frequency digrtal !'9"^^J.°J4!^ and processed by the receiving filter, 

obtained a complex base band signal I5.Q5 ''^"S ^^"^'^^^^^ is a signal obtained by subtracting the 

That is the signal by which the signal read out of he '^«"^° y'fJe,^Siate-frequency signal and correcting its initial 
larrfefirequen'cy error from the frequency of a '"^-"y.^^^tte Jt" to x ..x^.^o) 

phase, in other words, this signal is exp(-i(a-co , ''^^^J^^^^^^^ ^ime that is the sum of a deviation x^ from the 
fnstance, the initial phase is the sum of ^^P^^^/ "^^"^^^^ of the initial phase 60 by the carrier fre- 

symbol decision point and a deviation x^ ^^-^^^^^^^ Pf ^.^^^^^^^^^^ of the memory 1 3 to the base band signal has he 
cLncy error. The low-pa- fi«er for t^n^^^^^^^^^ ,e measured starting at the 

same charactenstic as that of the receiving tiiter a v example. 

beginning of the PN code has a length corresponding to 64 by 20 ^np orthogonal -transformed so as to 

^ T^en, an error vector is calculated from the ^^'^'^'^f^^^.^f^trca^t^^ modulation accuracy, amplrtude accu- 
corr^the frequency error arxi so forth and the ^^^^^^^^^^^^ quality p can be obtained 

racy phase error and IQ origin offset. S^ff'^'ent accuracy for me c^^^^^^^ 

Se base band signal 1,0^ but ^^^^f a'^^^^^^^ ^^"^'^ *J 

accuracy. Hence, the amplitude error, "^f ' Pj;Sio^3f,^ur?cy ^^^^ standard, and these estimated parame- 

rarsi«-=^^^^^^^^^ 

°nSrt:S-sepa.me,er.ar.es.^^^^^^^^ 

tude estimating section 24. Letting the number ^"^^^^'^"f^^^^f s^nal I5.Q5 by a(k) and e(k) and ttie 
ideal signal R,.Q, by I. the amplitude and Pj^!^ '^°"^°"^Q"^,;'eS aR(k), a carrier frequency error A<* the .nrtial 
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AcD = {(i;a(k) • [e(k) - e R(k)]) • (5:a(k) • k) - (Sa(k)) • {xa(k) • (e(k) - 
e R(k)] . k)}/{(5:a(k)) ' . (Ea(k)) • (^a(k) • k')} 

5 Co = {(5:a(k) . k) • (Ea(k) • [e(k) - e^{k)] - k) - (za(k) [e(k) - 

e R(k)]) • (za(k)}/{(5:a(k) • k)^ - (za(k)) • (5:a(k) • k^)} 

b 0 = {5:a(k) • a R(k) . cos[e R(k) - G(k) + Acd • k + G o])/5:a(k) ^ 

where £ indicates the summation from 1 of k to the last k (=N) of the input signal. 

These calculated values Aax Gq and bo are used to correct the Input complex base band signal Is.Qs by the following 
equation In a frequency, phase amplitude correcting section 25. 

75 

I 6 + jQ 6 = b o(l 5 + iQ 5)exp(-j{Aa> • k + 0 q) 

The thus corrected base band signal le.Qe and the Ideal signal R|.Q| are used to conduct the following equation in 
an error vector and waveform quality p calculating section 26. 
20 In the first place, squares of the amplitude of the error vector at symbol points (every chip period) are all added 

together by the following equation. 

r|c(k)|^ = 5:[(le(k). R,(k))' +(Qe(k). RQ(k))'] 

25 where T. indicates the summation from k=1 to 

This value Is divided by the number N of additions and the square root of the divided value is computed, followed 
by calculating the following equation to express it as a percentage with respect to the amplitude of the ideal signal. 

= 100 X V(X|£(k)|^/N)%rms 

30 

Next, to obtain the IQ origin offset, added values of the error vector in the I- and Q-axis directions are calculated by 
the following equations. 

5;£,(k)=£(l(k)-R,(k)) 

35 

XGQ(k)=£(l(k).RQ{k)) 

where Z indicates the summation from k=1 to N. 

The added value in each axis is divided by the number N of additions to obtain the mean value, which is converted 
40 to a dB value. Since the amplitude of the Ideal signal is 1 , letting the amplitude be represented by 0 dB, the IQ origin 
offset is given by the following equation. 

Oo = 20 X logio{[V{5:e,(k))%(SeQ(k))2]/N2} 

45 Next, the amplitude error is calculated. This begins with subtracting the amplitude 1 of the ideal signal from the 

amplitude of the measured signal le.Qe- followed by squaring the remainder. 

Zc^(k)2 = S(V{l6^(k)+Q6 2(k))-1)2 

50 where 52 indicates the summation from k=1 to N. This is converted by the following equation to the %rms unit to compute 
the amplitude en-or. 

£^,==100^x V(Ec;„2(k)/Ny- - 

55 After this, the phase error is calculated. In this instance, the phase of the ideal signal Ri.Rq is subtracted from the 
phase of the measured signal le.Qe and the remainder is squared. In this subtraction of phase, the phase difference 
may sometimes be calculated in the opposite direction when the phase is near the 180** axis; so that the phase differ- 
ence is corrected when it is larger than 180°. 
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Xu„ ^k) = s(arctan(Q(k)/l(k) - eR(k))' 

: converted by the following equation to a unit deg rms. 

2/ 



10 



15 



20 



25 



30 



35 



40 



45 



SO 



55 



P 

Where X indicates the summation from k=1 to N. This is ( 

Cp^ = (180Mx^(Xcp"(k)/N) 

The waveform quality p is calculated by the foilowing equations based on its definition. 

Re(i) = ^06M-RR(k) + Q6W-RQM) 
lmG) = S(Q6(k)-RRM-l6(k)-RQW) 

where x indicates the summation from k=l +640-1) to 64+64G-1). 

p = 5:(ReO) • ReG) + ImG) ' irnO)/[64 • X(l e + Q e \m 

obtain a time alignment error X3. . , On and T3 obtained as described above are d'spjayed on a 

The signal \s.Qe and the parameters p. J2L°°3o ^ade up of the orthogonal transform section 14 through 

be measured in the processing sections 23 *:°"f ?!. , ^ t^e input signal in the PN synchronized ideal signal 
While in the above ^ntersymbol interference oTdaS other than the symbol decision point may some- 

generating section 18. a signal with ^^^^'^^'^^^^^^^^Z not that of the CDMA digital cellular telephone^ In 
Les be needed in the case where held Ji a memory Although the present invention has 

such an instance. PN patterns processed by ^^^^/"f^ j''^^^^^^^^^^^^ quality of the modulated signal of the spread 
been described as being api^ied to the '^f!!" f "^""J ^veform quality of modulated signals of 

spectrum system, the invention is also applicable to «ord portion such as a synchronization 

^S^ersysteL-thiscanbeachi^edusu^^^^^ 

^ercSe%^r^"e^Sil^^^^^^ 

%SLalb«i above, according to the present invention. s^^^^^^^^ ^^t be 

^ ra^d7mru.r^^^^^^^^ — ^ 

input Signal, for example, about 400 symbols t^^^^^Jf^^^^^^ ideal reference signal thus obtained and 

erence signal is easily derived from the corrected '"P"*/'^^' '^^^^^^g^^) ^.^ors relatively easily. Thereafter, the input 
the correSed signal l^.Q^ are used to obtain .^^^"f^J. J^^^^^^^^^ signal k.O, so that these errors are 

mentioned U.S patent. Repeating the '=«'<="f'°"°*P?,''^2^e present invention, howe/er, the input signal .s sub- 

iSed-rrogo^^^^^^^^^ 

atun^ilation of calculation errors by repeating calculations. 
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By obtaining the frequency error, timing error and ideal reference signal with a signal shorter than the input signal 
to be measured, the computation complexity can be further reduced. 

In the case of performing the filter processing, the operation therefor is particularly time-consuming, but the com- 
putational complexity needed in the present invention can be reduced far smaller than in the prior art which always uses 
the entire signal to be measured. 

In the erhbodiment described above, since the frequency error is corrected prior to the receiving filter processing, 
part of the signal is not cut off by the receiving filter processing, ensuring accurate measurement. 

Since the CDMA mobile communication signal is spread by the PN code, if each chip is used as the symbol to be 
measured in the conventional method, the number of symbols becomes, for example, six-fold except the pilot signal and 
the computing time for measurement also becomes six-fold. However, the present invention is far smaller in computa- 
tion complexity as described previously, and hence it is advantageous over the prior art. 

According to the present invention, the input n-phase PSK signal is converted to a complex base band signal and 
it is raised to the n-th power as shown in Fig. 3. Accordingly, the carrier frequency of the complex base band signal is 
an error frequency of the input n-phase PSK signal, and this is a significantly small value and the sampling frequency 
of the signal raised to the n-th power satisfies the condition of the sampling theorem without fail, making it possible to 
obtain the carrier frequency with accuracy. 

In the synchronization of the spreading code with the input signal, the symbol decision point is estimated and the 
correlation is calculated for only samples near the symbol decision point as described previously with reference to Figs. 
4 and 5. the computational complexity required is smaller than in the case of calculating the correlation for every sam- 
ple. In the case of the pilot signal, in particular, since the correlation with the samples near the symbol decision points 
is calculated for only the real or imaginary part of the input complex base band signal, the computational complexity is 
sharply reduced as compared with that in the prior art. 

In the operation in the signal delay estimating section 22, that is, in the estimation of the symbol timing, the products 
A. B and C of the value of cross correlation between the input complex base band signal Z(k) and the ideal signal R(k) 
and the filter coefficients a^, b^ and are computed as indicated by Eqs. (16), (17) and (18) and t is calculated by 
Eq. (20). Accordingly, the value t can be obtained with high accuracy without involving any interpolating operation and 
hence in a short time (with small computational complexity). In contrast to this, according to the conventional method of 
Fig. 1, the calculation of the value t with high resolution and consequently with high accuracy, requires the reduction of 
the sampling interval by performing an interpolation operation. That is, the conventional method repeats the calculation 
while shifting the input signal and the ideal signal relative to each other, and hence it consumes much time for the cal- 
culation. 

It will be apparent that many modifications and variations may be effected without departing from the scope of the 
novel concepts of the present invention. 



1, A modulation parameter measuring method, in which: a digital sequence of an input digital quadrature- modulated 
signal stored in a memory is read out therefrom and subjected to orthogonal transformation into a complex base 
band signal; an ideal reference signal of said complex base band signal, synchronized with its symbols, is gener- 
ated; a frequency error, an Initial phase error and a timing error of said complex base band signal are estimated; 
said frequency error, said initial phase error and said timing error are used to correct said complex base band sig- 
nal; and said corrected complex base band signal and said ideal reference signal are used to measure modulation 
parameters of said input digital quadrature-modulated signal; 

CHARACTERIZED BY: 

A first step of reading out part of said digital sequence from said memory and subjecting it to orthogonal trans- 
formation into a first complex base band signal; 

a second step (30) of estimating a frequency error (hereinafter referred to as a first frequency error), an initial 
phase error (hereinafter referred to as a first initial phase error) and a timing error (hereinafter referred to as a 
first timing error) of said first complex base band signal and an ideal reference signal; and 
a third step (23) of obtaining said corrected complex base band signal by orthogonal transformation of the 
entire digital sequence necessary for measurement stored in said memory in such a manner as to correct said 
first frequency error, said first" initial phase error and said first timing error: ' ' ' 

2. The method of claim 1 . wherein said second step comprises: a seventh step (1 8) of generating said ideal reference 
signal by using said first complex base band signal; and 

an eighth step (19) of estimating a frequency error hereinafter referred to as a second frequency error) and 
Initial phase error (hereinafter referred to as a second initial phase error) from said first complex base band signal 
and said ideal reference signal, said second frequency error and said second initial phase error being respectively 
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said first frequency error and said first initial phase en-or in said third step. 

said first complex base band signal being said first initial phase error. 
The method of daim 2, wherein said seventh step compr^es: 3'^^°' J.^^^ ^^^^^^^ 

The method of claim 4. wherein said symbol decision point ^^^^^^^st^^^^^^^^^^^ 
lating the magnitude of each sample of said first complex base band s'Sna ■ a s^ep 0 35) °f Perf^^^^^ 

processing. 

The method Of daim 5. wherein said discrete Fourier transform processing is Pe^°-jf ^^y^-^'^P'^'"^ 
tuie sequence by instantaneous values of sine and cosine waveforms for each sample penod. 

. Themethodofdaim2.whereinsaKfsecondstepa>mprisesanin^^^^^^^^^^^^ 

eighth step being said first frequency error in said third step. 
The method of daim 7. wherein sad ninth step (15) --P^J-^ 

Signal of n-phase modulation to the n-th power: a step (1 P ^^^^^^^^ ^.equency com- 

:Se^^ireSrj^^^^^^^^ 

obtain said third frequency error. 
40 third step. 

beween said fct o«,pleK tese band sianal oorredrt .n ^^^'J^^^^J^^^l eh„acl»iste eaen 

^,;is;."is°rsr.sss:raS^^^^^^ 

results of filter processing using said three Nyquist filter charactenstics. 

X ■ o -J A «; R qand 10 wherein said second Step further comprises: a thirteen 
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between the address of said memory corresponding to said trigger and the address corresponding to the beginning 
of said ideal reference signal. 

14. The method of claim 14, wherein part of said digital sequence is about 400 symbols long. 

15. The method of any one of claims 1 through 14, wherein modulation data of said Input digital quadrature-modulated 
signal is preknown. 

16. The method of any one of claims 1 through 15. wherein said input digital quadrature-modulated signal is a trans- 
mitted signal modulated by a pilot PN code in a CDMA mobile communication system. 
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(54) Method for measuring modulation parameters of digital quadrature-modulated signal 



(57) The signal to be measured is converted into 
digital form and read into a memory 13 using a trigger 
as the reference, and a signal of a length 1/5 that of the 
signal to be measured, in the vicinity of the trigger, is 
orthogonal-transformed (14). A frequency error co^ of 
the transformed output li,Qi is estimated (15) and is 
used to correct the output l|,Q| (16). The corrected out- 
put \2,02 'S subjected to processing by a receiving filter 
(17), and the filtered output Is.Qs is synchronized with a 
FN code to obtain an ideal signal R|.Rq and a phase dif- 
ference Ti is calculated. The signal Is.Qa and the ideal 
signal Ri.Rq are used to obtain a frequency error 
and a phase error t2 of the signal la.Qs- All signals to be 
measured, stored in the memory 13, are subjected to 
orthogonal transformation (23) so that the frequency 
error to-i+cog and the phase error -c^+xg are removed, 
and the transformed output 15,05 and the signal R|.Rq 
are used to calculate modulation parameters as is the 
case with the prior art. 
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